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^ ' Abstract. We present a relativistic model for the motion of charged particles in 

in 

rotating magnetic field lines projected on to a plane perpendicular to the rotation 



axis. By making an approximation that the projected field lines are straight, an 



' analytical expression is obtained for the particle trajectory. The motive behind 

in '. 

^ developing this model is to elucidate some of the effects of rotation in pulsar 

fH ^ profiles. There is a significant contribution to the curvature of particle trajectory 

due to the rotation of pulsar, which is in addition to the inherent curvature of 
the field lines. The asymmetry in the observed pulse shapes can be explained by 



, considering the aberration-retardation effects. The single sign circular polarization 

that has been observed in many pulsars, might be due to the relative orientation 



of sight line with respect to the particle trajectory plane. 

Key words, pulsars-radiation mechanism: non-thermal radiation 



1. Introduction 

The wide diversity in the radiation characteristics of pulsars makes it difficult to under- 
stand fully the emission process in the light of models, which have been developed with 
some simplifying assumptions. Among the several emission mechanisms, curvature emis- 
sion has emerged as the most probable emission mechanism (Sturrock 1971; Rudcrman 
& Sutherland 1975, hereafter RS75; Lyne & Manchester 1988; Gil & Snakowski 1990). In 
order to explain the high brightness temperature observed in pulsars, coherent emission 
by bunched of particles has been postulated (Karpman et al. 1975; RS75; Buschauer & 
Benford 1977). Also, other models based on plasma effects have been proposed for pulsar 
radiation (e.g., Mehose & Gedalin 1999; Asseo & Rozele 2000; Gil et al. 2004). 

Most of these models give emphasis on explaining the high brightness temperature of 
pulsars, but the polarization is poorly explained. However, the polarization observations 
such as the polarization angle swing favors the curvature radiation. It has been considered 
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as a natural emission process for pulsars, though there are unresolved problems like 
the bunch formation, orthogonal polarization modes etc (e.g., Stinebring et al. 1984; 
Gangadhara 1997; Gil et al. 2004). 

It is imperative to understand the influence of rotation to study closely the curvature 
emission mechanism. The idealized case of particle acceleration has been discussed by 
Machabeli & Rogava (1994, hereafter MR94), who considered particles moving freely 
along an infinitely long, rigidly rotating straight tube and derived an expression for the 
trajectory of a particle. Gangadhara & Lesch (1997) have proposed a model for the 
particle acceleration in rotating magnetosphere in the context of active galactic nuclei 
(AGN). Reiger & Mannhiem (2000) have also discussed the particle acceleration along 
the rotating straight magnetic field lines in AGN, by assuming the angular velocity of 
particles is same as that of AGN. 

In the case of pulsars. Gold (1969) was the first one to propose a pulsar emission mech- 
anism based on rotation. This model was further taken up by many authors and found 
to encounter difficulties in explaining the interpulses (e.g., Sturrock 1971). Blaskiewicz 
et al. (1991) have studied the effects of corotation velocity on the pulsar radio emission by 
assuming a constant emission height. Hibschhman & Arons (2000) have extended their 
work to include the first order effects to study the delays in the phase of polarization 
angle sweep due to aberration. Later, Peyman & Gangadhara (2002) have improvised the 
model of Blaskiewicz et al. (1991) by relaxing the assumption of constant emission height, 
and analyzed the effect of rotation on the morphology of pulsar profiles and polarization. 

Gangadhara (1996, hereafter G96) has derived the equation of motion of a charged 
particle in pulsar magnetosphere. He has considered the straight field lines, which are 
projected on to a two dimensional (2D) plane placed perpendicular to the rotation axis. 
The dominant forces, which act on a particle moving along rotating the field lines, are 
the magnetic Lorentz force, centrifugal force and coriolis force. The rotational energy of 
the pulsar is transferred to the corotating plasma as it moves along the field lines. The 
magnetic Lorentz force acts as a constraining force and drags the plasma along the field 
lines. Because of the inclination of magnetic axis relative to the rotation axis, corotating 
plasma tends to rotate with an angular velocity less than that of pulsar on some field 
lines. The difference in the angular velocities of particle and pulsar had been already 
pointed out in RS75. 

In the present work, as a follow up of G96, we consider the same 2D geometry and 
analyze the dynamics of a charged particle. Since the field line curvature radii of open 
field lines are comparable to the light cylinder radius, over a significant radial distance, 
we can approximate them to be straight lines. 

In super strong magnetic fields, the drift velocity becomes negligible compared to the 
velocity parallel to the field lines. The Larmour radius of gyration becomes quite small, 
and hence particles almost stay on the same field lines all along their trajectories. This 
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motion is considered as the bead-on-wire approximation. The particles are accelerated 
because of the unbalanced centrifugal force, and thus extract the rotational energy of the 
pulsar. The single particle emission is considered in this model and plan to consider the 
collective effects in our follow up works. We take into consideration of the non- uniform 
angular velocity of particles, which can be less than the pulsar angular velocity on some 
field lines which are inclined with respect to the meridional plane. Since the particle 
trajectories are found to be curved, we estimate the curvature emission and analyze the 
effects of rotation on the radiation characteristics. In § El and § 13 we solve the equa- 
tion of motion of a relativistic charged particle and find its trajectory. We compute the 
characteristic frequency of curvature radiation in §0] In §|31we estimate the polarization 
parameters and plot them with respect to different parameters. 



2. Charged particle dynamics 

We assume that the dipolar magnetic field lines are projected onto a plane perpendicular 
to the rotation axis. Consider an inertial Cartesian coordinate system as shown in Fig.^ 
where the 'z' axis is parallel to the rotation axis ($7) of pulsar. The projected magnetic 
axis on the x-y plane coincides with the x-axis at time t = 0. The equation of motion for 
a charged particle moving along a rotating magnetic field line is given by (G96), 

m— I = m\l r, (1) 



dt\ dt^ 

where m — mo7 is the relativistic mass, 7 the Lorentz factor, toq the rest mass, $7* the 
angular velocity and r the radial position of a particle. 

Let Vr = dr/dt and = rQ* be the components of particle velocity, then 

13 =^{Vrer + ¥4,64,), (2) 

where c is the speed of light. We define the unit vectors in the radial and azimuthal 
directions as 

Br = {coscj), sincj), 0) , (3) 
60 = (-sine/), cos(/), 0) , (4) 

where is the angle between the radial vector to particle and the a;-axis. Then the 
Lorentz factor of particle is given by 

-1/2 

(5) 



1 

7 



IdrV frVL* 



c dt ) V c 



Consider a particle injected at the point B onto a magnetic field line which is inclined 
by an angle 0p with respect to the a;-axis at time t = 0. Let do ~ OB be the distance 
between B and the rotation axis. The effective angular velocity (G96) of a particle is 
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Fig. 1. The coordinate system in which the particle motion is considered. The curve BQ 
represents the particle trajectory in the x-y plane. 



given by 



— (Iq cos^ 9o — do sin Oq \/ — dt^ cos^ (?o 
r — do' cos^ ^ c^o sin 



62 



(6) 



where f2 is the angular velocity of pulsar, h = (iocos6'o and 6*0 = (7r/2) — 0p is the angle 
between the field line tangent and at B. Using the relation for f2*, we can write 7 as 



7 



1 + 



1 dr 
c dt 



-1/2 



where I? = r2c?o cos^q/c Thus, using the expression for 7, we rewrite Eq. JQ): 



d2r d7 dr 2 m 



7 r 



By multiplying Eq. ((SJ by rji^c^), and defining a dimensionless variable 
_Vt r 

we obtain 



fs [2s^~Dyil + D'')] fds 



dt^ 



\dt 



1 + £|2 



(7) 



(8) 



(9) 



(10) 
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Since 6*0 is close to 7r/2 for the field lines, which are close to the x-axis, we find <C 
for dp < ^- Therefore, we reduce Eq. ({Tn|l by dropping the terms containing + D^, 

and obtain 



To find the solution of Eq. Ijf 1(1 . we follow the method proposed by Zwillinger (1989). By 
choosing / = [ds/dt)'^, we can be reduce it to the following form: 

f + /^.2.tf. (12, 

Its solution is given by 

f ^n^{l-s^) + C{l-s^f, (13) 

where C is the integration constant. To find C we use the initial condition that at t=0, 
it follows from Eq. (j^Jl that sq — s(0) — rofl/{c\/l + D"^) and sq = ds/dt\t=Q = 



vt)Q,/{c\/l + -D^), where tq and vq are the initial position and velocity of particle. 
Therefore, we obtain C ~ — fi^fc^ and 
1 



1 - s2 
1 6q 



1- 



(l-.g)Si=. 

Hence Eq. (|13|) reduces to the following form 
ds 



(14) 



— = f}^(l - s2) - A:2(l - .s2)2 (15) 
whose solution is given by 

s = cn(A-r2t), (16) 
where cn(z) is the Jacobian Elliptical cosine function (Abramowicz & Stegun 1972), and 

sn(z) 

/ , 17) 

i ^(1-W2)(l_fc2y;2) 



/ , (18) 

^ 



(^) 

Using the expression for s given by Eq. © , we find the radial position of the particle 



= arccos f j . (19) 



cJl + 1)2 

r = cn(A - m) . (20) 

The radial position of the particle according to Eq. (|20() as a function of time is plotted 
in Fig. 121 It shows that the particle position increases with time and reaches a maximum 
at the distance of light cylinder radius = Pc/27r, where P is the pulsar period. Next, 
the particle returns back to origin, due to the reversal of the centrifugal force. This type 
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Fig. 2. Radial position of the particle as a function of time. Used 70 = 100, do = 10^ cm 
and 6*0 = 90°. 

of oscillatory motion of a particle in an infinitely long, straight and rigidly rotating tube 
has been discussed by MR94. 

Though we have extended the calculation of r of a single particle all the way up 
to light cylinder, it may not be realistic in the case of plasma motion. Near the light 
cylinder, plasma inertia causes the field lines to sweep back and break down of the rigid 
body motion. The Lorentz factor of a particle, which follows from Eqs. lO, 1)15(1 and 116|l . 
is given by 

7= ^ . (21) 

k^{l + D^)sn^{X-nt) 

3. Particle trajectory and its radius of curvature 

In Fig. Owe consider a particle moving along the field line BQ. The point A represents 
the particle injection point at time t ~ that is at a distance do from the rotation axis. 
The particle co-ordinates can be defined as 

{x, y) = r{t) (cos(/>tot, sin^tot) , (22) 

where (/<tot is the angle between the radial vector to the particle and the a;-axis. Based 
on Fig. |31 we define 



(23) 



where 




(24) 



For do <C r, we find 



(25) 




OA X 

Fig. 3. The geometry of motion of a particle along a rotating field line BQ. The angles 
are ZXOM = nt, ZMBQ = (j>p and ZMOQ = 0', and the radius OA=OB=rfo- 



The ± signs in Eq. (|23() corresponds to the sign of the angle (pp . In Fig. 0] we have 
plotted the trajectories of the particles moving along different magnetic field lines, which 
are marked with (pp. It shows that the trajectories are curved towards the direction of 
rotation of pulsar. The particles moving in those trajectories are accelerated, and hence 
they emit curvature radiation. The curvature radii of those trajectories slightly differ 
from one another, as the particle angular fl* is different for each field line. 

To derive the curvature radii of the particle trajectory, we approximate cn(A — Qt) 
and r{t) using the formalism given by Pearson (1974): 



cn{z, k) = cos z + —{z — sin z cos z) sin 2; + 0{k'^). 



(26) 



In the limit of i 1 and A: ^ 1, the series expansion of r(t) is given by 

r{t) = + ait + a2t^ + ai,t^ + a^t'^ (27) 



where ag, oi, 02, 03, 04... are the expansion coefficients: 

fc^ sin A (A — cos A sin A) " 



ao 



ai 



a2 



n 



cos A + 



16 



32 



2 sin A{8 - fc^ + 3 fc^ cos(2 A)} - 4 /c^ A cos A , 
16-9 k^) cos A + fc^ (9 cos(3 A) + 4 A sin A) 



03 = — — — ffc2 {4 A cos A + 13 sin A- 27 sin(3 A)} -16 sin A] 

96 



04 = ^ 7 [(16 -nP) COS A + fc2 {81 cos(3 A) + 4 A sin A}] 

384 



(28) 
(29) 
(30) 

(31) 
(32) 



For vo « c, Eq. lO implies fc w 0. Therefore, using Eqs. (HHJ, CZI), CHIl and l(T<^ . 

we find A = 7r/2 and sn z = sin z. Thus, we have 



r{t) w -^-^ sin (m) . 



(33) 



8 



Thomas & Gangadhara: Radio Emission by Particles due to Pulsar Spin 




0.02 0.04 0.06 0.08 0.1 0.12 
X /Tt 



Fig. 4. Particle trajectories during the time interval <t < 0.02 sec in laboratory frame. 
We considered the field lines which lie in the range of —15° < (pp < 15° with an interval 
of 5°. Assumed neutron star radius i?n ~ 10 Km. 



Using the expression for r{t) and Eqs. H22|l and H25|l . we find the curvature radius of 
particle trajectory: 

_ [{dx/dt)^ + [dy/dtff/'^ 

^ ~ {dx/dt){d^y/dt^) - {dy/dt){d^x/dt^) 

~ ^tlVi + D^. (34) 

It shows that the curvature radius of a particle trajectory is approximately /2 for 
do ~ 7r/2. However, for other values of we find p becomes slightly larger than rj^jl. 
Note that these values of p are comparable to the intrinsic curvature radii of dipolar field 
lines in the emission region given by Gangadhara (2004). We find, in the conal emission 
regions, the intrinsic curvature radii of field lines are comparable to the curvature radius 
induced by rotation into the particle trajectory. If the core emission is believed to come 
from the field lines which are close to the magnetic axis then it becomes difficult explain 
the core emission due to the intrinsic curvature of field lines. This is because the field 
lines which are very close to the magnetic axis have very large curvature radii, and for 
the magnetic axis it is infinity. Therefore, in the absence of rotation, we can not expect 
any significant curvature emission from the particles or plasma moving along the field 
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9 




lines which are close to the magnetic axis. However, from the observations we do see 
many pulsars emitting strong cores. In our model we show that if we consider rotation 
of particles, the core emission can be explained. Because the rotation induces significant 
curvature into the trajectory of particles which move along the field lines that are close 
to the magnetic axis. 



4. Radiation emitted by a particle 

When the particles move along curved trajectories, they emit curvature radiation. The 
characteristic frequency of the curvature radiation is given by (RS75) 

« 317 ^-^2. (35) 

By having known the 7 and p from Eqs. iQ and H34|) . we can estimate lOc- The total 
power emitted by a particle is given by 

3c fc4[l + i:)2-(r/rL)2]4' ^•^^> 

where q is the particle charge. In Fig. [S] we have plotted W as a. function of 0p at 
r = 0.5 r^. Wc find maximum power is emitted by the particles which move along the 
field line with 0p = as = 0, compared to those on the other field lines {(j)p ^ 0). 
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5. Polarization of Radiation 

The radiation electric field is given by (e.g., Jackson 1972; Gangadhara 1997) 

E(tj) = C/ J hx{nx /3)exp{iw(t-h.r/c)}dt, (37) 

— oo 

where C/ = —iuj q e^'^^°^'^ / V27rSoc, So is the distance from the origin to observer, uj 
the radiation frequency and n the sight line. To solve the integral, we shall express 
h X in X (3) (see, Appendix-A) and the argument of exponential as series expansions in 
time t. Consider the sight line which makes an angle 9 with the 2D plane, and 77 with 
the x— axis: 

n — (cos cos 77, cos sin 77, sin 6*). (38) 

To describe the polarization state of the emitted radiation, we define orthogonal unit 
vectors (see. Fig. P): 

e|| = (— sin^cosTy, — sin6'sin?7, COS0) , (39) 

e_L = (— sinr/, cos?7, 0) . (40) 

The unit vectors (ri, e||, e^) form an orthogonal triad: 

n X = e|| . (41) 

Let Iq be the time at which (3 aligns with fi, and the observer receives radiation. We 
transform the time variable t to ^ + such that fJto stands for an initial phase. Thus, 
we find 

r{t + to) = a'„ + a[t + a'^t^ + a'^t^ + a'^f^ ... (42) 

based on Eqs. and l(?7jl . The expansion coefficients Og, a'^, 03, 04 ... are same as 
ao, ai, 02, ... except A replaced by A — ^IIq. 

Using Eqs. (|20|l and (|22|l . we find the series expansion of the exponential argument 
in Eq. H37|l and keep the terms up to the order of t^: 



(t + io) 

c 



T 

t + to COS 9 (cos 77 cos ^tot + sin 77 sin (f>tot) 

c 



No + Nit + N2t^ + N^f .... (43) 



The coefficients No, Ni, N2 and N3 are given in the Appcndix-C. The series expansion 
of exponential argument is converging, and it is quite obvious from Eqs. 120(1 and ((26|l . In 
the limit of fc w 0, the series expansion of r behaves like the trigonometric sine function. 
Since the angular width of emission beam is ~ 2/7, the time taken by the particle to 
cross the angular width of the order of emission beam is sa 2p/c'y. Thus the truncation 
of higher order terms introduce a negligible error in our calculations. Since we intend to 
reduce the integral in Eq. (|37|l to a known form, we limit the series expansion terms up 
to the order of t^. 
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Using the transformation given by Buschauer & Benford (1976), we find the electric 
field components (see, Appendix-B): 

£|| = ^ (V^ioBo + V^|i-Bi + VpB2) Cf e^^« sinO , (44) 

E± = \ iV±oBo + V^iBi + Vi_2B2) Cf e'^° . (45) 

Next, we define the Stokes parameters as 

I = + E_lE*^ , 

Q = -B||£'|* — E_\_E\ , 
[/ = 25R , 

V = 2'=s [e^Es^ . (46) 
The linear polarization is given by 

L=^/Q^TTP. (47) 

In observations, pulsar polarization is normally expressed in terms of L and V. So, in 
our model, we call L and V as polarization parameters. 

5.1. Polarization parameters of radiation emitted by many particles 

We consider a set of field lines on the 2D plane, and estimate the total emission by par- 
ticles accelerated along them. During pulsar rotation, the sight line stays at a particular 
with respect to the 2D plane. Since the emission from each particle is relativistically 
beamed in the direction of velocity /3, the observer tends to receive the radiation from 
all those particles, for which (3 falls with in the angular width of ±1/7 with respect to h. 

First we estimate the polarization parameters of the radiation emitted by a single 
particle at the instant to < fmax. The instant to is the time at which n.l3 = 1 for a given 
initial ^p. As the rotation progresses, new to is computed for the advanced rotation phase 
by again solving n./3 = 1, and computed the polarization parameters. This procedure is 
continued till to ~ ^max, where tmax is the time at which the particle goes out of radio 
emission zone {do < r < 3 x 10^ Km). Since the radiation is emitted over a range of r, 
due to the aberration and retardation the radiation beam gets shifted to the leading side 
of the pulse. The role of retardation and aberration phase shifts has been discussed by 
e.g., Phillips (1992) and Gangadhara & Gupta (2001). 

In order to compute the total polarization parameters with respect to the rotation 
phase, first we sort the polarization parameters due to single particles into groups of 
phase bins and add them. In the following steps, we give the details of the procedure: 

1. Fixed the observer's sight line at a specific 6 with respect to the 2D plane. 

2. Selected a set of field lines in the range of —5° < < 5° with a successive line 
spacing of 0.1°. 
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3. Solved n.p = 1 to find to Sit the point of emission on the trajectory corresponding to 
each field line, and estimated the polarization parameters at those points. 

4. The retardation phase shift ri(imax — to) is subtracted from rj assigned for each of the 
emission beam, and estimated the effective rotation phase. 

5. Next, the sight hne is rotated by 0.1° to a new phase, and repeated the procedure 
(1-4) over the range of -12° < 77 < 12°. 

6. Finally, the array of polarization parameters are sorted into groups of phase bins, and 
added to get the total pulse profile. 

In Fig. |B| we have given the total polarization parameters computed from the emis- 
sions by many particles as functions of rotation phase. In panel (a) we have plotted the 
profile that is obtained when the sight line lies in the 2D plane, and panel (b) for the case 
when the sight line is inclined by —0.05°. The profiles indicate that the peak emissions 
are shifted to the earlier phase as a consequence of the aberration- retardation effect. 

6. Discussion 

For simplicity we considered the dipole field lines projected on to a 2D plane. The rel- 
ativistic particles are assumed to be streaming outward along such a field configuration 
with an initial Lorentz factor in the range of 10^ < 70 < 10"^. Since our aim is to under- 
stand the rotation effects on particle dynamics and pulse profile, we consider the single 
particle emission, and leave the collective plasma emissions to latter works. We are inter- 
ested in the region which extends from a few stellar radii to a radial distance well within 
the light cylinder, where the radio emission is expected to occur and the bead-on-wire 
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approximation holds. Though we have approximated the particle motion to 2D, we are 
able to investigate the influence of rotation on pulsar profiles. Our model is more relevant 
for the cases where the inclination angle a of the magnetic axis relative to rotation axis is 
large enough. In such cases, the projected field lines may be approximated to be straight 
lines over a significant radial distance. We have derived an expression for the radial po- 
sition of a particle fEo. I2U|) . which shows an oscillatory behavior, as shown in Fig.|21 A 
similar case of particle motion in an infinitely long, straight and rigidly rotating tube, 
has been discussed by considering a gedanken experiment by MR94. They have shown 
that due to the centrifugal force reversal, particle returns back to the rotation axis after 
reaching a maximum distance at which the rotation velocity reaches the speed of light. 
It turns out to be an oscillatory motion in the radial direction. Gangadhara (1996) has 
shown that the particle angular velocity cannot be same as the field line angular velocity 
if the magnetic axis is inclined with respect to the rotation axis. We considered this 
effect in our treatment of particle motion, and found the particle trajectories and their 
curvature radii vary with field line orientation. 

Since the magnetic field is very strong, Larmour radius of gyration and drift velocity 
becomes very small. So, the particles are assumed to follow the same set of field lines 
all along their trajectories. In the case of single particle dynamics, the magnetic field 
dominates and hence the rigid body motion may be extended all the way up to the light 
cylinder. But in reality plasma corotates with the neutron star, and we must to take into 
account of the plasma inertia in the region close to the light cylinder. Therefore, it is 
possible that the magnetic field lines will sweep back, and can lead to the generation of 
toroidal magnetic field. Hence the oscillatory motion that our solution predicts, can not 
be achieved in a real physical situation like pulsars. So, the particle which reaches the 
vicinity of the light cylinder can not come back, but escape from the magnetosphere as 
a pulsar wind. 

We find the energy of particle increases due to the centrifugal force, as indicated by 
the Eq. for Lorentz factor 7. By this way the rotational energy of neutron star gets 
transferred to the particles via the magnetic field lines. 

We find the radius of curvature of particle trajectory is approximately ri/2, which is 
comparable to the inherent radius of curvature of dipolar field lines (Gangadhara 2004). 
So, we believe the curvature emission due to the rotational motion of particles should be 
comparable to the actual curvature emission in corotating frame. Both the Lorentz factor 
and the characteristic frequency reach the maxima at the region close to light cylinder. 
Hence our model indicates that the high frequency radiation (e.g.. X-ray, 7-ray) may be 
emitted in the regions close to light cylinder. 

In a later work, Rogava et al. (2003) have shown that if a particle freely moves along 
a tube with an arbitrary curvature, the centrifugal force does not reverse always. They 
have showed that the particles move in the tube with a variable angular velocity. This 
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supports our result that the particles angular velocity on some field lines differs from 
that of pulsar. That is the particles moving along the field lines with (j)p — rotate with 
the angular velocity which is same as the pulsar angular velocity. But those moving along 
other field lines, for which cjjp ^ 0, rotate with the angular velocity which is smaller than 
the pulsar angular velocity. The particles moving along the field lines with (f>p ^ 0, tend 
to emit more power than those moving along other field lines, and hence the profile in 
Fig. [5l shows a peak at <j>p ^ 0. Also, it is evident from observations that the peak of 
pulsar profiles (core) is, probably, emitted from the field lines with (f>p ^ 0. 

By taking into account of aberration-retardation, we have reproduced a simulated a 
pulse profile (Fig. O by adding the radiation emitted by particles accelerated on a set 
of field lines. The sign of (j) ba-s been flipped to match with the phase sign convention 
followed in pulsar profiles. The roughness in the curves of Fig.Elare due to the increments 
of 0.1° in (j)p and Tj. This choice was made based on the limitation in computing time. 
However, the smoother profiles can always be generated by choosing smaller increments 
and opting for longer computing time. Since we consider an uniform plasma flow along 
the field lines, our profiles do not have subpulsc components. 

Our model shows the effects such as the aberration and retardation makes the 
pulse profiles to become asymmetric about the pulse center. This phenomenon has 
been observed in most of the pulsar profiles (e.g., Gangadhara & Gupta 2001; Gupta 
& Gangadhara 2003). 

In our model, we find if the sight line is at a fixed angle {6 = —0.05°, see. Fig. [SJ to 
the particle trajectory plane, observer tends to receive a single sign circular polarization. 
This type of single sign circular polarization has been observed in many pulsars (e.g., 
Han et al. 1998). As a foUowup of this work, we plan to consider the full 3D dynamics 
of plasma in a rotating dipolar magnetic field, and estimate the coherent radiation. 



7. Conclusion 

By considering projected dipolar magnetic field lines on a plane perpendicular to the 
rotation axis, we have developed a 2D model for the particle dynamics in pulsar magnc- 
tosphere. The motive behind developing this model is to elucidate some of the rotational 
effects induced in the pulsar profiles. We have obtained the analytical expressions for 
the particle trajectory and its curvature radius. The energy of particles increase at the 
expense of neutron star's rotational energy. We find the sight line orientation relative 
to the particle trajectory plane might determine the sign of circular polarization. The 
asymmetries observed in the pulse profiles can be explained by considering the aberration- 
retardation effects. 
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Appendix-A: To find a series expansion for the factor n x (n x /3) that appears 
in Eq. (|T7|> 

Consider f3 in Cartesian co-ordinates: 

/3=hv.x + Vyy), (A-1) 

where x and y are the unit vectors along the x and y-axes, respectively, (Fig. 0. Then 
it follows from Eq. (O that 

= V, COS I j -V^ sm I 1 , (A-2) 

and 

Vy = V, sm ( j COS ( 1 . (A-3) 

Using Eq. (|42|l . we derive the series expansions for radial velocity and rotation velocity 

= rQ* . (A-5) 
By substituting Vr and into by Eqs. ljA-21) and ljA-3|l . we obtain 

= VyO + Vylt + Vy2t\.. . (A-6) 

The expressions of V^q, Vyo, V^i, Vyi ... in the above expansions are lengthy (see, Mathew 

and Gangadhara 2005). Using the triple vector identity and the definitions of n, e|[ and 
e^, we obtain 

n X (n X /3) = -(/3 . e,,) e,, - (/3 . e^) , (A-7) 
where 

sin 9 , , , . ^ 

l3.i\\ = ^ (Fy sin?/ + Ux COST?), (A-8) 

f3.e± ^ -{Vy cos -V^ sin ri) . (A-9) 
Using the series expansions of and Vy, we write 

hx {nx ^) = ^ [e\\sme{V\\o + V\\it + V\\2t^) + i±iV^o + V^it + V^2)] , (A-10) 
where 

V\\ i = Vyi sin ry + T^i cos rj , (A-11) 

Kli = V^i sin 77 — COS 7/ , (A-12) 

and i = 0, 1, 2 . 
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Appendix-B: Transformations for solving Eq. (|37|> 

Using the method of Buschauer & Benford (1976), we make the following transformations 
in order to solve the integral in Eq. 1)37(1 . 
Consider 



1 



exj)[i{Nit + N2t^ + N^t^)] dt = — TTre*^" / exp[i(zr + r^)] dr 



where 



is a dimensionless variable, and 
1 



N. 



1/3 



Nl 



By differentiating the Eq. I)B-1|I with respect to A^i and N2, we obtain 

00 

1 



t exp[i{Nit + + N^r)] dt 



N. 



2/3 



J T exp[z(zr + r^)] dr 

—00 
00 

Ci / exp[i{zT + T^)] dr 



and 



00 

j t'^ exp[i{Nit + 



+ iVat^)] dt 



N. 



1/3- 



Cm J exp[z(zr + t^)] dr 

— OG 

00 

Cp / r exp[i(zT + r^)] dr 



We define 



Li(z) = y exp[i(zT + r3)]dr-^V^Ki/3[2(z/3f/2], 



-00 
00 



exp[z( ZT + t3)] dr = z-= z K2/3[2(z/3)3/2] 
V 27 



Bo ^ 

Bl : 
B2 : 

where 

a = 

Cn = 



1 


e^^ 






1 


e*^ 






1 


e'^ 







*^"[L2(z)-CiLi(z)] 



^e'^"[GnLi(z)-CpL2(z)] , 



N2 



3N.^ 
N2 
37V3 



2/3 ' 



2Nl 
97V3 



iVi 



(B-1) 



(B-2) 



(B-3) 



(B-4) 



(B-5) 

(B-6) 

(B-7) 

(B-8) 
(B-9) 
(B-10) 

(B-11) 
(B-12) 
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Cm — 



Co 



-ml 

9Ni 
2N. 



3N3 



4/3 
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(B-13) 
(B-14) 



Appendix-C: The expressions for expansion coefficients which appear in 
Eqs. (in 10 ) 

The expression for (fitot in Eq . H23|) is expanded in a series given by 

(/-tot =00 + 01^ + 02 i^ + 03t^ 



t>o — arccos/i3 , 



Ail Oq 



io sin'' 0p ^c?o'' | cos(4 — cos(2 <^p)^ 
+ flo' do^ |l + 5 cos(2 (/)p)| flo' + 2 ao''^ — 2 do /Ui |rfo^ cos(3 
+ 3 COS0P ao'^l ) { ^ 2 ai'^ + 4 ao' a2'| 



Mio 1 ^^i/ii3cos0p do^i2sin^0p^ 



3 v^l - V 3 ao' 
The coefficients No,Ni, N2 and which appear in the series expansion of Eq. (|43|l 
are given by 

cos 9 cos(?7 — 0o) ao' ' 



Ni — — [c — cos 9 {(1)1 sm{rj — 0o) clq' + cos{ri — cjio) oi'}] , 
c 

N2 = ^ cos 6* {01^ cos(?7 - 0o) - 2 02 sin(77 - (/)o)} ao' 
2 c L 

- 2 sin(77 - 0o) ai' + cos(77 - (/)o) a2'} 



iVs = — cos 6* 
6 c 



|6(?!)i 02 cos(ry - 0o) + (0i^ - 6 03) sin(77 - 0o)| ao' 
+ 3 |0i^ cos(77 - 0o) ~ 2 02 sin(77 - 0o)) ja/ 
- 6 |0i sin(?7 - 0o) a2' + cos(?7 - 0o) as'j . 

The terms appearing in the series expansion in ljA-6|l are given by 



/2 ,2 

ao - 



1/2 



Vxo = - yfl sin0oao' \J oq'^ - 6^ ) +cos0oai 

Vxi = —01 sin 00 ai' + 2 cos 00 02' 

+ 01 cos 00 (^b^ — ao'^^ ~ sin 0o ao' ai 

K2 = — [01^ cos 00 + 2 02 sin 0o] a/ — 4 0i sin 0o 02' + 6 cos 0o 03' 

+ 2 cos 00 {b — ao') (6 + ao') {(ao'^ — b^ )02 + +0i ao' ai'} 



+ Q sin 00 {b* 0i^ — 2 6^ 0i"^ ao'"' + 0i^ ao'^ + 



2 J, 2 /2 



2 „ /4 



(C-1) 
(C-2) 

(C-3) 



(C-4) 

(C-5) 
(C-6) 
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V^o = ^ cos \J ao'^ — IP' + sin ai' 
Vyi = (pi cos (f)o ai' + 2 sin 00 ^2' 



+ 



bi Ct sin (j)o (b^ — ao'^^ + O cos (po do' ai' (^^o'^ — b 



-1/2 



sm 00 



y2 = 3 sin 00 as' + (02 cos0o ^ 

+ { ~ ("Ai^ cos 00 + 2 02 sin 00^ (b"^ - ao'^) 

+ 2 01 sin 00 oo' (b^ — oo'^^ ai' — b^ cos 0o ai'^ 

+ 2 cos 00 ao' (ao'' - 6") 02'}] [2(ao'' - b^^^ 



a\ + 2 01 cos 00 02' 
2 



Ml = A/l- 



j 2 • 2 
rfo sin 



ao 



/2 



/"2 
/"3 



do sin 



■y^l - (/X2 + /Ui cos 0p)^ 



/2 2 

/U4 = ao /ii 



Ms = 
Me = 

M7 = ( 

Ms 
M9 



/U2 cos 0p aj^ 



\i\ Oq sm 0p 

/ mIqi y 

Vjuf ao sin^ 0p/ 



M2 



/xi sin0p 
/xi /X7 cos 0p 



3 a? 




/2 
"0 


^ - Me 


^a'l^ 




.a? 





Mio 



Mil 



M12 



Ml3 



= {(I + 2/Z32) /Z52-2/Z3 (m3^-1) Ms 

+ |c?o (1 + 2 /X3^) sin^ 0p ai' (/U2 a'l -2/^5 ao') | ^ | - — ^ — -3 , 

' a(i' > (/i32 - 1) 

(-4 /is /X8 + (-1 + Ms^) M9) (m5 ao'^ - rfo sin^ 0p ai'^ 

2(l-M32)^ao'' 

-6 (^2 ai'^ - 3 ao' ai' 02' + ao'^ a/^ 



ao 



/3 



ai'^ - 2 ao' ai' 02' + ao'^ as' 



ao 



/3 



M2^ (Mii Qq' + 3/f7ai') 
/if 2 ao' sin^ 0p 
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